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The final state — black hole or naked singularity — of the gravitational collapse of a marginally 
bound matter configuration in the presence of tangential stresses is classified, in full generality, in 
terms of the initial data and equation of state. If the tangential pressure is sufficiently strong, config- 
' urations that would otherwise evolve to a spacelike singularity, result in a locally naked singularity, 

' both in the homogeneous and in the general, inhomogeneous density case. 
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; I. INTRODUCTION 

T— I . 

^ ! The final state of gravitational collapse remains as one of the outstanding problems in classical general relativity. 
For a given set of Cauchy data, hyperbolic evolution of the field equations leads to a unique final state and, under 
certain conditions, a singularity is formed However, it is not known whether this singularity is spacelike, or a 
visible spacetime singularity. Understanding how such singularities arise from regular initial data, and whether they 
—i ' can be visible (at least to local observers), is still to a large extent an open issue 

—i , The most studied analytical solution is that describing spherical dust collapse, whose detailed analyses over the past 
^ ' ten years have provided many valuable insights into the formation, visibility, and causal structure of singularities 
^ Albeit of considerable interest in its own right, dust collapse constitutes a highly idealized model, both in terms 
of geometry and matter content. Whilst the former may arguably be a reasonable approximation for astrophysical 
^ collapse [Q, the existence of pressures in high density regimes, together with an effective equation of state, cannot 
bJ[)> be neglected in realistic situations Q. It would, therefore, be highly desirable to understand the final state of 
' spherical gravitational collapse for generic equations of state. Unfortunately, this objective looks still quite far from 
being reachable, in particular in the case of perfect fluid sources. Therefore, research has turned to less ambitious 
objectives. An obvious way to generalize dust (i.e. stress-free) models, is to try to "add" non-vanishing stresses. 

A model in which analytical treatment appears feasible is that of vanishing radial stresses, since in this case the 
general exact solution is known in closed form Q (the opposite case, in which only a radial stress is present, has 
also been recently considered |8|). The formation and nature of singularities in gravitational collapse with tangential 
stresses has recently been studied quite extensively [|[ |l^, In the models analyzed so far, a tendency to 

uncover part of the singularity spectrum was observed, i.e., initial configurations that would otherwise end up in black 
holes develop locally naked singularities (in one special case, the singularity was even shown to be timelike |^, in 
contrast with the naked singularities of inhomogeneous dust collapse, which are typically null). However, all such 
models belong to a very special subclass of solutions with tangential stresses, the so-called Einstein cluster. This is a 
system of counter-rotating particles, in which the stress is generated by angular momentum. From the physical point 
of view, the Einstein cluster is interesting since it mimics the effects of rotation without introducing deviations from 
spherical symmetry. However, although such a system can be formally obtained by the choice of a particular function 
of state within the general exact solution, this function is essentially a "centrifugal potential" and therefore does not 
fulfill the physical characteristics which are typical of the state equations of matter continua, like those to be expected 
in strongly collapsed matter states (e.g. in neutron stars). As a consequence, the results obtained for the Einstein 
cluster cannot be straightforwardly extended to the general case of tangential stresses, which is of course interesting 
in its own right. 
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Thus motivated, in the present paper we carry out a complete analysis of the final state of gravitational collapse 
with tangential stresses for (marginally bound) realistic matter configurations. We show that the final state depends 
on the first non- vanishing term of the Taylor expansion of the initial density distribution near the center (as happens 
to be the case in dust models) and on the first non- vanishing term of the Taylor expansion of another function which 
carries all the relevant information coming from the choice of the material. It follows that the final fate of the collapse 
depends on the choice of two integers, n and k, giving the order of these two expansions, respectively. As we shall 
see below, there is a critical value for k, above which the tangential pressure effects are negligible, and the endstate 
of collapse is indistinguishable from that of dust models. Below the critical value, however, tangential stresses come 
into play, and all the configurations that would otherwise end up in a black hole, terminate in a singularity that is 
at least locally visible. At the critical parameter, a transition behavior is observed, wherein the singularity may be 
visible, depending on the details of the initial data. Both homogeneous and inhomogeneous density evolutions are 
examined in detail, and they are found to be qualitatively equivalent. 

This paper is organized as follows. In Sec. II, we present a short but self-consistent account to recall the physical 
framework and to prepare the needed mathematical formalism. The results are given in Sec. Ill, where the structure 
of the spectrum of endstates is presented and discussed. The paper ends with concluding remarks in Sec. IV. 

Geometrized units, in which G = c = 1, are used throughout. 



II. SPHERICAL COLLAPSE WITH TANGENTIAL STRESSES 



A. The general solution 



In this section we review briefly the main properties of spherically symmetric gravitating systems with vanishing 
radial stresses (for details, see 0]). 

The mathematical structure of the field equations is simplified by the fact that the Misner-Sharp mass is conserved. 
This actually allows complete integration in the so called mass-area coordinates, in which the line element reads 



ds^ = -K^ (1 ]dm^ + 2 — 

R J uh 



dRdm ~ ^dR^ + R^dn^ 



(2.1) 



Here K, u and h are functions of R and m, and dVP' = dO^ + sin^ OdcjP' is the canonical metric of the unit two-sphere. 
The function u is the modulus of the velocity of the collapsing shells and satisfies 



2m 

~R^^~h? 



while 



where 



K = g{m) + J G{m, R)dR , 



1+2 U 



2 m 



/l2 



-3/2 



(2.2) 



(2.3) 



(2.4) 



The functions g{m) and E(m) are arbitrary. The quantity h{m, R) is a measure of the internal energy stored in the 
material, and can thus be regarded as an effective equation of state (the physical characterization of this function will 
be given below). It can be shown that the energy density e and the tangential stress 11, can be written in terms of h 
as: 



h 



n = -— 



AttuKER'^ ' 
R dh 

2h'dR^' 



(2.5) 
(2.6) 



These formulae show that the tangential stress vanishes whenever the function h is independent of R. In such a case, 
the material is not sustained by any internal stress and the line element ( p.l| ) reduces to the one describing spherically 
symmetric dust in mass-area coordinates [|l^. Since, in all formulae, only the ratio of E and h appears, the value 
of h along an arbitrary curve R = Ro{m) can be rescaled to unity, so that, in particular, the dust solutions can be 
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characterized hy h = 1. The function Ro{m) plays the same role as that played by the initial mass distribution in 
the standard comoving coordinates [the inverse transformation from mass label to comoving label being r — RQ^m)], 
while the function E{m) is the specific binding energy function. 

The above recalled structure shows that a solution with tangential stresses is identified — modulo gauge 
transformations — by a triplet of functions {g,E,h}. This parameterization contains the dust spacetimes as the 
subset {g, E, 1}. In this way, we can construct in a mathematically precise fashion, a comparison between dust and 
tangential stress evolutions, by comparing the end state of the dust collapse {{g, E, 1} with chosen g and E) with the 
endstates of the tangential stress solutions {g, E, h} with the same g and E, and different choices of the equation of 
state h. 



The physical singularities of the spacetime described by the metric (2.1) correspond to infinite energy density, and 
are given by i? = or if = 0. At i? = 0, the shells of matter collapse to zero proper area, thereby leading to 
shell-focusing singularities, whereas the vanishing of K implies intersection of different shells of matter, corresponding 
to shell-crossing singularities. If occurring, such singularities are gravitationally weak and we shall not deal with them 
in the present paper. 

Non-central shells which become singular (R — and m ^ 0) are always spacelike — thus covered — in spherical 
spacetimes with radial stresses, as a consequence of mass conservation |Q; only the central singularity R = m = 
can be visible, since at R — m ^ 0, the apparent horizon and the singularity form simultaneously. The visibility is 
determined by the existence (or lack thereof) of outgoing null geodesies with past endpoint at the singularity. To 
investigate the existence of such geodesies, we use the general method originally developed by Joshi and Dwivedi p[ . 
One defines the quantity x — i?/2m" (where 1/3 < a < 1) and observes that, at the singularity, both R and m vanish. 
Using I'Hopital's rule, one can arrive to an equation which, in the present model, takes the form 



xo — lim 



m2 



R->0 



2a 



Ro,mh{m, Rq) 
E{m)u{'m, Rq) 



GdR 



Ra 



E^ 1 
(-1 + tt)™""^ + - 

/l^ x 




E^ 



(2.7) 



If this equation admits real positive-definite roots, i.e. at least one finite, positive solution xq exists for some a, then 
this root represents the tangent to an outgoing null geodesic meeting the singularity in the past, which is therefore 
at least locally naked (covered otherwise). 



B. Equation of state, regularity, and physical reasonability conditions 

As recalled in the previous section, the choice of the matter model in a solution with tangential stresses corresponds 
to the choice of the function h. This is equivalent to the choice of an equation of state like, e.g., the barotropic 
equation of state connecting density and pressure in the case of a perfect fiuid. The choice of h is thus restricted by 
physical considerations, as follows [ p^ . 

The internal energy per unit volume of the material (e, say) is given by e = pbh, where pb is the baryon number 
density. If the material is physically viable, the internal energy has to coincide with the number density in a (locally) 
relaxed state of the continuum. This obviously means, as in the classical (i.e. non-relativistic) mechanics of continua, 
that h must be a positive, convex function having a minimum in the locally relaxed state. Without loss of generality, 
we take the value of the minimum equal to one. Since all deformations are described "gravitationally" in the frame 
we are working with, the state of local relaxation corresponds to the flat space values of the metric components. It 
follows that the most general physically valid equation of state has the form 

h{m, R) = l + (3{m){Ro - Rf + /i2(m, R), (2.8) 

where h2{m, R) goes to zero more rapidly than {Rq — R)^, as R tends to Rq. The function /3(m) is positive and models 
the "strength" of the tangential pressure. With the scaling Ro{ra) = r, the squared factor on the right-hand-side lies 
in the interval [0, r^], where the lower and upper limits are realized on the initial and singular slices, respectively. 
Since we are interested only in the behavior near the central singularity, it follows that, as far as the causal structure 
of the singularity is concerned, we can consider the function h to be of the form 

/i(m, R)^l + Pkfn^'^iRQ ~ Rf , (2.9) 

where k must be positive in order to assure regularity at the center, and the positive constant (5k measures the effects 
of tangential stresses. The case of dust collapse corresponds to = 0, while for (3k non-zero the first non vanishing 
term of the Taylor expansion of the tangential stress near the center is of order k. 

Having chosen a suitable matter model, to address the issue of cosmic censorship, one must also ensure that several 
other reasonability conditions are satisfied. 
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First of all, the spacetime must have a regular center (so that the singularity forming at r = is a true dynamical 
singularity). This requires the functions if, u, and E/h to be smooth and bounded away from zero. Further, we 
require the collapse to proceed from a regular initial slice, where the metric and the exterior curvature are continuous, 
the physical quantities are bounded, and there are no trapped surfaces. As we have seen, initial data consists of 
two independent quantities g and E, which can be equivalently described by the initial mass [m{r), i.e., RQ{m)] and 
the initial velocity distribution / := E^ — 1. In what follows, we take f — 0, i.e., we consider marginally bound 
configurations. Regularity on the initial surface implies that 

m(r) = For^ + F„r''+^ + C'(r"+'^), (2.10) 

where Fn denotes (47r/3 times) the first non- vanishing derivative of the initial energy density profile at the origin. 
Positivity of mass requires Fq > 0, and Fn < for any realistic configuration, where the energy density decreases 
away from the center. 

We want to exclude the presence of trapped surfaces on the initial slice. This requires 2m < i? on the Ro{m) = r 
surface: 

2m\ 2m(r) „ , „ ,o 



R 



= 2For-' + 0{r''+^) <1, (2.11) 

R=r ' 



where the second equality holds only near the origin. The first equality imposes a maximum size on the matter 
distribution, for a given density profile. Near the origin, the requirement 2m/ R < 1 is always satisfied on the initial 
slice, for < r < 1/V2K- 

Finally, and in accord with the spirit of the cosmic censorship conjecture p^ , we require that the matter content 
obeys the weak energy condition (WEC): 

e > , e + n > 0. (2.12) 



From Eq. (2.5), one sees that the first inequality is sat isfied provided h > and K < (recall that m < 0, for the 
collapsing situation we are interested in). From Eq. (|2.9| ), it follows that dji{\nh) < 0, and hence the tangential stress 
is positive, which is a sufficient condition for the second inequality (positivity of 11 implies that also the strong energy 
condition, e + 211 > 0, is automatically satisfied here). 

C. The final state of collapse 

Dust models do not admit globally regular solutions, i.e., singularity-free solutions. This is not necessarily true for 
general matter fields, even within spherical symmetry. Examples of globally regular solutions in spherical symmetry 
include certain classes of perfect fluid models, which evolve to an eternally oscillating configuration from regular 



initial data (see e.g. |17|), and solitonic solutions of the massive Einstein-Klein- Gordon system [|18|. Accordingly, it 
is important to understand the qualitative behavior of the dynamics, from a given set of regular initial data, before 
addressing the causal natur e of the final state. 



The equation of motion (2^) for the collapsing shells can be re- written as 



2 ^ Z{R, flo) 

R[l + P{m){Ro - R)^] ' ^ ' ' 

where the "effective potential" Z is defined by 

Z{R, Ro) := 2m + f3{m){2m - R){Ro - Rf . (2.14) 

This function may be regarded as the analogue of the Newtonian effective potential governing the motion of the fixed 
shell i?o, wherein the allowed regions of the motion correspond to Z > 0. Setting R — ({t, Ro)Ro, to the lowest order 
in Rq we obtain 

2Fo[l + hm^'^R^ - Rf] > hm'^'^it). (2.15) 

This condition always holds, independently of the choice of k, as to goes to zero. Therefore, there is always a region 
of initial data leading to continued gravitational collapse. 

We now would like to check whether the central point (m — 0) eventually gets trapped. Since m is conserved for 
any given shell, this can be done by analyzing the dynamics of shells near the center, along curves R ~ Xm with 
A > 2 j|l^. The requirement Z >0 now reads 

2> Phm'^^^iRQ - Xmf{X-2) . (2.16) 



5 



Clearly, near the center this inequality always holds, which implies that non-central shells inevitably become singular. 
This is in contrast with the Einstein cluster case, where rotation has the effect that all the shells near the central one 
remain regular for all times [|lOj. In the present case, instead, the apparent horizon inevitably forms, and therefore 
the central singularity may or may not be visible, depending on the existence of geodesies meeting this singularity in 
the past. 



III. THE SPECTRUM OF ENDSTATES 



A. The root equation 



The analysis of the root equation for gravitational collapse with tangential stresses is non-trivial, because of the 
integral appearing on the right-hand-side of Eq. ( [3.5[ ), which introduces a "non- local" behavior. Indeed, as we have 
seen, only the case of the Einstein cluster has been analyzed so far. 

To check whether the root equation has positive-definite solutions, we observe that, although the integral appearing 
in square brackets cannot be carried out in terms of elementary functions, its behavior near the center can be analyzed 
as follows. First, let us denote 



h{m;a) 



2m" a; 



G{m,R)dR, 



(3.1) 



-Ro 



where G{m,R) is defined by Eq. (2.4). Due to the mean value theorem, there exists x("*) ^ {Ro,2m°'x) such that 

/i(m; a) := (2m"x - Ra)G{m, x(m)). (3.2) 

Since Rq ^ m^l'^ as m goes to zero, 1/3 < a < 1, and x is positive-definite, it follows that both to~^/'^x("^) ^-i^d 
vnf^ l \(vfi) are also positive-definite as m 0. One can thus evaluate the right hand side of Eq. ( |3.2D — using Eq. 
(2.4) — as follows: 



/i(m;a) = 



2\/2 



1/2 



X //5[F~^/^-(x/mi/3)]2^ 



ii/3 I 2[l + /3(i?o-x)^ 



-3/2 



1 



X /(m/3,™)[i^-i/^-(x/mi/3)]\ ^ (/3F-i/^)[^^-i^=^ - (x/mV3)] 



ml/3 \^ 2[1 + /3(i?o - x)']' 
[l + /?(i?o-x)]'/', 



,1/3 



[1 + /?(i?o - x)' 



(3.3) 



where the dots stand for terms of higher order in positive powers of rri^l^ . Since /?(0) — 0, the divergence is driven 
by 0(mr^\ with all the other terms remaining finite. Therefore, the limit 



/2(?n;a) :— lim / mG{m,R)dR, 



(3.4) 



is convergent. Using Lebesgue's dominated convergence theorem, we can expand the integrand near the center (m = 0) 
in leading powers of m, and integrate successive terms | p^ . It follows that the root equation, for a general density 
distribution [given by Eq. (2.10)], with the equation of state (2.£), can be written as: 



™|(1-q) / 9™3(a-l)/2 



^0 



a 



where 



1 TO 

Jx 



(3.5) 



e„ := 



Bk 



nFn 



(2n-|-9)/6 ■ 



18a/2Fo^ 
\/2(9 + k)(3k 



945iv 



3/2 



(3.6) 
(3.7) 
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B. The dust limit 



In order to better visualize the effects of tangential stresses in spherical dust collapse, it is important to recall briefly 
how to recover the (very well known) dust case from our general set-up of the previous section (for a complete review 
and list of references see [^). This is easily done by setting /?(m) = 0, and thus (3k = 0. The root equation becomes 
then 



xq = lim — 

m^o a 



e„m(3+2n-9a)/6 ^ 



,3/2 



+ . 



-1/2 



1 ~ X ' m 



l/2^(l-a)/2 



(3.8) 



It is straightforward to check that for a 7^ a„ = i(l + ^) the root equation does not admit any positive-definite 
solution: if a < q;„, then = 0, and for a > the right-hand-side of Eq. ( |3.8|) diverges in the m — > limit. One 
must therefore have a = a„. In this case, the equation reduces to 



3/2 



3/2' 



1 - x-^'^ lim m(3-«)/9 



(3.9) 



Clearly, we must have n < 3, for xq G (0, -f 00). For n < 3, we obtain 



Xq 



2n 



2/3 



(3.10) 



and the singularity is locally visible. For n = 3, after a little algebra, one obtains the quartic equation (where 
Z = y^): 



with 



2^-* + Z3 + 7(1 - Z) = 0, 



= 03 = /{2V2) 



(3.11) 



Standard results from polynomial theory can be used to show that this quartic has real and positive roots iff 7 > 7c, 
where [|j 



:= (26-H 15V3)/4. 



(3.12) 



C. Structure of the endstates 

For (3k ^ 0, the root equation is given by Eq. (^^). It is immediately seen that one can always choose 1 < n < 3 
and/or 1 < ^ < 3, in such a way that, for suitable 1/3 < a < 1, this equation gives a finite, algebraic condition. It 
then follows that in all such cases the singularity can be naked. 

For n < 3 or k < 3 the singularity is naked, since a positive root always exists (we omit here tedious, but 
straightforward calculations; the values of the roots are reported in Table I). 

For fc > 4, the effects of tangential stresses are negligible, and the endstate is exactly the same as that in dust 
collapse, with the shape of the initial central density profile determining the existence of positive-definite roots, and 
thus local visibility. The very same can be said, however, about the effects of inhomogeneity for n > 4, since all such 
configurations would lead to black holes, while now they all terminate in a visible singularity if A; < 3. 

A transitional ("critical") behavior is observed at the three entries forming the "boundary" of the n < 3, k < 3 
region. Here, naked singularities occur iff a quartic equation has positive solution(s). This quartic is identical in form 
to that holding for dust, namely 

2Z^ + + 7„fe(l - Z) = 0, (3.13) 

where 7„fc takes the values 

733 = 7 + -^, (3.14) 

734 = 7, (3-15) 
743 = S, (3.16) 
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with 



105Fo 



3/2 ■ 



(3.17) 



Therefore, at each "boundary" , naked singularities are formed iff the corresponding quantity ^nk is greater than the 
quantity 7c defined in (3.12). 



Initial data n=l n—2 n=3 n > 4 



k=l 
k=2 



-21Fi+8/3if|P^ ^13^ ^13^ ^13^ 

SAV2fI^"' 1^ 1^ 1^ 

-Fi -105F2+22feFp^/^ life life 

4,V2fI^'^ 420V2fI^'^ 210V2fI'^ 210V2F^^^ 

T^F^ I7#7e 733 > 7c 743 > 7c 

k>4 — ""^^i /R — ~^?^/K 734 > 7c no positive real roots 

~ 4,V2Fg '° 4,V2Fg'° 

TABLE I: The endstate of collapse for inhomogeneous collapse with tangential pressures. See text for details. 



D. The homogeneous case 



The homogeneous dust solution is the well-known Oppenheimer-Snyder spacetime ||2C| (the seminal paradigm for 
black hole formation, for over sixty years), whose dynamical evolution results in a Schwarzschild black hole, with an 
inter ior, spacelike (hence covered) singularity. To evaluate the root equation for this case, we simply set = in 
Eq. (3.8), obtaining: 



1/2 



— lim — 

m-*o 3a 



(3.18) 



Clearly, we need a < 1, or else the right-hand-side of the above equation diverges negatively in the m ^ limit. For 
a = 1, we obtain Xq^'^{2xq^^ + 1) — 0, which has no positive-definite solution, and the singularity is therefore covered. 
For a < 1, the Xq terms drop out, and a self-consistent solution exists iff a = 1/3, which is not allowed. Hence, the 
singularity is always covered, as expected. 

This example provides a particularly simple test-bed fo r st udying the effects of tangential pressure on the final state 
of collapse. In presence of tangential stresses, from Eq. (3.5) with 8„ = 0, we obtain: 



lim — 



,(l-")/2 



(3.19) 



A curious phenomenon occurs here. In fact, the above equation is formally identical to the root equation for inho- 
mogeneous dust [cf. Eq. (3^) above] with k and Bk playing the role of n and On respectively. Clearly, in the leading 
term we can always choose 1 < fc < 3, such that 1/3 < a < 1, and the singularity may thus be visible. Again, there 
is a clear transition at fc = 3, below which (stronger tangential stresses) the singularity is always visible, and above 
which (weak er ta ngential stresses) is always covered. At the transitional value, fc = 3, we have again to solve a quartic 
of the kind ( 3.13| ) with parameter equal to S, so that naked singularities occur if (S > 7c. The final state of collapse is 
summarized on Table II, below. 



Initial data Root (a;o 



,3/2 



Singularity 



k = 1 
fc = 2 
fc = 3 



21F, 

11^2 



3/2 



5 > 7c 



visible 
visible 
transition 



fe > 4 no real positive roots covered 



TABLE II: The endstate of homogeneous collapse with tangential pressure. 
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IV. CONCLUDING REMARKS 



We have examined here the effects of tangential stresses on the final state of (marginally bound) collapse. Since 
the nature of the singularity turns out to depend only on the Taylor expansion of the data and of the state function 
near the center, it is possible to classify the final state of all the solutions that have a physically reasonable equation 
of state. This classification can be done in terms of the "strength" of the tangential stress near the center. If this 
strength is above a certain threshold, the tangential pressure dominated collapse always leads to a locally naked 
singularity, independently of the choice of the initial data. If, instead, the tangential pressure is sufficiently weak, 
then the final state is the same as that occurring in dust collapse, where the shape of the initial density profile near 
the origin fully determines the visibility of the singularity. 

These results actually disproof a conjecture, recently put forward by one of us Q, the proposal of which, roughly 
speaking, was that the space of the data leading to naked singularities in spherical collapse could not be "larger" 
than that of dust initial data leading to the same singularities. Contrary to this expectation, the present results are 
consistent with quite a different scenario: they support the notion that the strong version of the cosmic censorship 
conjecture is likely to be violated for realistic matter distributions in a way that can safely be described as common, 
insofar as the matter content is concerned. Indeed, our analysis shows that naked singularities are not an artifact 
of dust models, or a byproduct of "toy" models constructed without giving the equation of state of the collapsing 
matter, or, finally, a spurious outcome of unnatural effects induced by rotation in spherical symmetry. 

There remains the possibility that such naked singularities are an artifact of spherical symmetry by itself, although 
both numerical and analytical results suggest that dynamical curvature singularities may persist in non-spherical 
collapse H, ||. 
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